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Abstract 

The exclusive production of p° mesons in deep inelastic electron-proton scattering has 
been studied using the ZEUS detector. Cross sections have been measured in the range 
7 < Q 2 < 25 GeV 2 for 7*2? centre of mass (cm.) energies from 40 to 130 GeV. The 
7 P — > P P cross section exhibits a Q 1 ' ' -°- 5; dependence and both longitudinally and 
transversely polarised p s are observed. The ^y*p — > p p cross section rises strongly with 
increasing cm. energy, when compared with NMC data at lower energy, which cannot be 
explained by production through soft pomeron exchange. The data are compared with 
perturbative QCD calculations where the rise in the cross section reflects the increase in 
the gluon density at low x. the gluon density at low x. 
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1 Introduction 



With the high energy electron-proton collider HERA, it has become possible to study deep 
inelastic scattering (DIS) processes at large values of Q 2 , the negative of the four-momentum 
transfer squared of the exchanged virtual photon, and large values of W, the virtual photon- 
proton centre of mass (cm.) energy. The exclusive production of vector mesons in DIS is of 
particular interest. While numerous data exist at low Q 2 [1-5] on the exclusive reaction 

e (or jj) + p — > e (or p) + p° + p, (1) 

only two experiments have reported DIS measurements for Q 2 > 5 GeV 2 |6|, |7|. These latter 
measurements have been restricted to W < 20 GeV. 

Previous studies of exclusive leptoproduction (7*iV — > p°N) and real photoproduction 
(<yN — > p°N) off a nucleon N have shown that for p° production at low Q 2 (typically < 2 
GeV 2 ): 

• the Q 2 dependence of the cross section can be described by the Vector Dominance Model 
(VDM) [|TJ in which it is assumed that the photon fluctuates into a p° meson yielding: 



da(Q 2 
dF 




(2) 



where £ is the ratio of the longitudinal to transverse forward amplitudes, e is the relative 
longitudinal polarisation of the virtual photon and M p is the p° mass; the distribution of 
t, the square of the four-momentum transfer between the photon and the p°, is described 
by a single exponential dependence, in the range from t = to t — —0.5 GeV 2 , with a 
slope parameter, b ~ 7 — 12 GeV -2 ; 

in real photoproduction (Q 2 = 0), the process is 'quasi-elastic' and the helicity of the 
p° is similar to that of the incident photon, i.e. s-channel helicity is largely conserved 
(SCHC) [§]. The p° decay distribution exhibits an approximate sin 2 9h dependence, where 
9h is the polar angle of the tt + in the p° cm. system and the quantisation axis is the p° 
direction in the 7p cm. system; 

the real photoproduction p° total cross section increases slowly as a function of W for 
W > 15 GeV H |TD[. This is expected for an elastic reaction dominated by the exchange 
of a 'soft' pomeron with an intercept of the Regge trajectory of a(0) = 1 + e' = 1.08. The 



intercept is determined from fits JTTJ] to hadron-hadron total cross sections: for a 7p total 
cross section a tot ~ ^ 2 (a(o)-i) _ W 2e'^ the optica i theorem yields ^ _ n ~ W ie ' ~ W - 32 . 



t=o 



At larger Q 2 (2 < Q 2 < 25 GeV 2 ), leptoproduction results from the EMC @ and NMC 0] 
experiments indicate that: 

• the 7*p — > p°p cross section is consistent with a l/Q 4 behaviour; 

• at larger Q 2 (> 6 GeV 2 ), the distribution of the square of the transverse momentum of 
the p° with respect to the virtual photon (p\) is exponentially falling with a slope of 
b = 4.6 ± 0.8 GeV~ 2 [0], about a factor of two smaller than that of the photoproduction 
elastic process; 
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as Q 2 increases, the fraction of zero-helicity, longitudinally polarised p°s increases above 
50%; assuming SCHC [p]], this implies that the longitudinal virtual photon cross section 
dominates; 



for 2 < W < 4 GeV the cross section falls with W at small Q 2 < 4 GeV 2 
significant dependence on W is observed for 9 < W < 19 GeV M. 



No 



The reaction ^*p — > p°p has also been the focus of theoretical investigations. Early studies 
based on VDM are described elsewhere jl| . A study of diffract ive leptoproduction by Donnachie 
and Landshoff based on a zeroth order perturbative QCD (pQCD) calculation for pomeron ex- 
change at small values of Bjorken x ~ Q 2 /W 2 [[11]] reproduced many of the features seen 
experimentally, including the Q 2 dependence of the data at low W. In more recent presenta- 
tions, the pomeron is treated as a non-perturbative two-gluon exchange fL2"| . This approach 
has also been studied by Cudell ||13|| . Calculations in pQCD have been performed in the leading 
logarithm approximation for J/ip electroproduction by Ryskin [14]. Ginzburg et al. [15] and 
Nemchik et al. |16| have also performed a calculation in pQCD for vector meson production. 
These more recent calculations predict a Q~ 6 dependence of the longitudinal cross section at 
high Q 2 , in contrast to the Q~ A VDM expectation. Brodsky et al. |17 have studied the forward 



scattering cross section for this reaction by applying pQCD in the double leading logarithm 
approximation (DLLA). They predict that at high Q 2 the vector mesons should be produced 
dominantly by longitudinally polarised virtual photons with a dependence for the longitudinal 
cross section: 



da L 



dt 



t=o 



A 
Q~ 6 



l + i(7r/2)( 



d 



d In x' 



xg{x,Q 2 ) 



(3) 



where A is a constant, which can be calculated, and xg(x, Q 2 ) is the momentum density of the 
gluon in the proton. Using the x-dependence of xg(x, Q 2 ) measured at HERA and the relation 
W 2 ~ Q 2 /x for small x, at fixed Q 2 one expects that ^ ~ W 1A [fT8| , in contrast to the 

W ' 32 dependence expected for the 'soft' pomeron. 

This letter presents a measurement with the ZEUS detector of the exclusive cross section for 
p° mesons produced at large Q 2 by the virtual photoproduction process 7*p — > p°p at HERA. 
The data come from neutral current, deep inelastic electron-proton scattering in the Q 2 range 
of 7 - 25 GeV 2 , similar to that of the earlier fixed target experiments || [7|; however, they cover 
a lower x region (4 • 10 -4 < x < 1 • 10 -2 ) and, consequently, a higher W region (40-130 GeV). 



2 Experimental conditions 

The experiment was performed at the electron-proton collider HERA using the ZEUS detector. 
During 1993 HERA operated with bunches of electrons of energy E e = 26.7 GeV colliding with 
bunches of protons of energy E p = 820 GeV, with a time interval between bunch crossings of 
96 ns. For this data-taking period 84 bunches were filled for each beam (paired bunches) and 
in addition 10 electron and 6 proton bunches were left unpaired for background studies. The 
electron and proton beam currents were typically 10 mA. The ep cm. energy is y'i = 296 GeV 
and the integrated luminosity was 0.55 pb _1 . 

ZEUS is a multipurpose magnetic detector whose 1993 configuration has been described 
elsewhere [fL9|. This brief description concentrates on those parts of the detector relevant to 



the present analysis. 
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Charged particles are tracked by the inner tracking detectors which operate in a magnetic 
field of 1.43 T. Immediately surrounding the beampipe is the vertex detector (VXD) |2(J which 
consists of 120 radial cells, each with 12 sense wires. Surrounding the VXD is the central 
tracking detector (CTD) |2l] which consists of 72 cylindrical drift chamber layers, organised 
into 9 'superlayers'. In events with charged particle tracks, using the combined data from both 
chambers, resolutions of 0.4 cm in Z and 0.1 cm in radius in the XY planeQ are obtained for 
the primary vertex reconstruction. These detectors provide a momentum resolution given by 
a PT /p T = ^J(0.005p T ) 2 + (0.016) 2 (with p T in GeV). 

The superconducting solenoid is surrounded by a high resolution uranium/scintillator calor- 
imeter which is divided into three parts: forward (FCAL), barrel (BCAL), and rear (RCAL) 
covering the angular region 2.2° < 9 < 176.5°, where 9 = 0° is defined as the proton beam 
direction. Holes of 20 x 20 cm 2 in the centre of FCAL and RCAL are required to accommodate 
the HERA beam pipe. The calorimeter parts are subdivided into towers which in turn are sub- 
divided longitudinally into electromagnetic (EMC) and hadronic (HAC) sections. The sections 
are subdivided into cells, each of which is viewed by two photomultiplier tubes which provide 
the energy and the time of the energy deposit with a resolution of better than 1 ns. An addi- 
tional hadron-electron separator (RHES)|[22||, located at the electromagnetic shower maximum 
in the RCAL and consisting of a layer of 3 x 3 cm 2 silicon diodes, was used to provide more 
accurate position information for electrons scattered at low angles than was available from the 
calorimeter alone. 

The luminosity is measured from the rate observed in the luminosity photon detector of hard 
bremsstrahlung photons from the Bethe-Heitler process ep — *> epj. The luminosity detector 
consists of photon and electron lead-scintillator calorimeters |p3|| . Bremsstrahlung photons 
emerging from the electron-proton interaction point at angles below 0.5 mrad with respect to 
the electron beam axis hit the photon calorimeter placed 107 m along the electron beam line. 
Electrons emitted at scattering angles less than 5 mrad and with energies 0.2E e < E' e < 0.9E e 
are deflected by beam magnets and hit the electron calorimeter placed 35 m from the interaction 
point. 

The data were collected with a three-level-trigger. The first-level-trigger (FLT) for DIS 
events required a logical OR of three conditions on sums of energy in the EMC calorimeter 
cells. Details are given elsewhere [|l^, |24|]. For events with the scattered electron detected in 
the calorimeter, the FLT was essentially independent of the DIS hadronic final state. The FLT 
acceptance was greater than 97% for Q 2 > 7 GeV 2 . The second-level-trigger used information 
from a subset of detector components to reject proton beam-gas events, thereby reducing the 
FLT DIS triggers by an order of magnitude, but without loss of DIS events. 

The third-level-trigger (TLT) had available the full event information on which to apply 
physics-based filters. The TLT applied stricter cuts on the event times and also rejected beam- 
halo muons and cosmic ray muons. Events remaining after the above veto cuts were selected 
for output by the TLT if 5 = Eji£j(l — cos^) > 20 GeV — 2E^, where Ei, 9i are the energy 
and polar angle (with respect to the nominal beam interaction point) of the geometric centre 
of a calorimeter cell and E 1 is the energy measured in the photon calorimeter of the luminosity 
monitor. For fully contained events 5 ~ 2E e = 53.4 GeV. For events from photoproduction, 
the scattered electrons remain in the rear beam pipe and S peaks at low values. 



1 The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of HERA. 
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3 Kinematics of exclusive /r production 



The kinematic variables used to describe p° production in the reaction: 

e + p -»• e + p° + X, (4) 

where X represents either a proton or a diffractively dissociated proton remnant of mass Mx, 
are the following: the negative of the squared four-momentum transfer carried by the virtual 
photon^ Q 2 = —q 2 = — (k — k') 2 , where k (k f ) is the four-momentum of the incident (scattered) 
electron; the Bjorken variable x = -£pr, where P is the four- momentum of the incident proton; 
the variable which describes the energy transfer to the hadronic final state y = |-^; the cm. 
energy, y/s, of the ep system, where s = (k + P) 2 ; W, the cm. energy of the 7*p system: 
W 2 = (q + P) 2 = Q + M 2 ps ys, where M p is the proton mass; and t' = \t — t min \, where 
t is the four- momentum transfer squared, t = (q — v) 2 = (P — P') 2 , from the photon to the 
p° (with four- momentum v), t m i n is the minimum kinematically allowed value of t and P' is 
the four-momentum of the outgoing proton. The squared transverse momentum p\ of the p° 
with respect to the photon direction is a good approximation to t' since t' is, in general, small 
(<< 1 GeV 2 ). For the present dataQ t min ranges from -0.0006 to -0.08 GeV 2 . 

In this analysis, the p° was observed in the decay p° — > 7r + 7r _ . The momentum vector 
of the p° was reconstructed from the pion momentum vectors determined with the tracking 
system. The production angles (9 P and 4> p ) and momentum (p p ) of the p° and the angles of the 
scattered electron (9 J and <p e '), as determined with RCAL and RHES, were used to reconstruct 
the kinematic variables x,Q 2 , etc. The energy of the scattered electron was determined from 
the relation: 

E c = (s + M 2 n - M 2 x )/2 - (E e + E P )(E P - \p p \cose p ) 
e ~ (E e + E p ) (1 - Pcos9' e ) - (E p - \p p \cos9 ep ) ' 1 ' 

where E p is the energy of the nir pair, 9 ep is the angle between the tttt three vector and the 
scattered electron and j3 = (E p — E e ) / (E p + E e ). For the case of reaction (1), Mx = M p and El 
is a good estimator of the energy of the scattered electron, E' e ; for events in which the proton 
diffractively dissociates into the system X, Mx > M p and El is only slightly different from E' e . 
The above expression simplifies to 

El ps [2E e - [E p - \p p \cos9 p )}/(l - cos9' e ) (6) 

when Mx = M p and the transverse momentum of the proton is negligible compared to its 
longitudinal component. This last relation provides an accurate way to calculate the kinematic 
variables and is a simple expression used to evaluate the radiative corrections for this process. 
The variable y is calculated from the expression y = (E p — \p p \cos9 p )/2E e . The calculation of 
p% also uses the p° and electron momenta: p\ = (p ex + p px ) 2 + (p ey +Ppy) 2 - 



4 Monte Carlo simulations 



The reaction ep — > ep°p was modelled using two different Monte Carlo generators. The first, 
DIPSI [p^l , describes elastic p° production in terms of pomeron exchange with the pomeron 



treated as a colourless two-gluon system \T%\. The model assumes that the exchanged photon 



2 In the Q 2 range covered by this data sample, effects due to Z° exchange can be neglected. 
3 The -0.08 GeV 2 value corresponds to M x = 8 GeV and Q 2 = 25 GeV 2 . 
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fluctuates into a quark-antiquark pair which then interacts with the two-gluon system. The 
cross section is proportional to the square of the gluon momentum density in the proton. 
Samples of and u events were generated in a similar way. 

A second sample of p° events was generated with a Q~ 6 dependence for the ep reaction, a 
flat helicity angular distribution, and an exponentially falling p\ distribution with a slope of 
6 = 5 GeV -2 . The Monte Carlo generator used the HERWIG framework [26] and the events 
were weighted according to the measured helicity, p\ and y distributions. 



A third p° Monte Carlo generator (RHODI), based on the model of Forshaw and Ryskin |27 



was used to model the proton dissociative processes with a da(j*p) / dM\ oc l/M^ 5 dependence. 
Different Mx dependences were obtained by weighting the events. Events were generated for 
M\ values between 1.2 and 4000 GeV 2 . All Monte Carlo events were passed through the 
standard ZEUS detector and trigger simulation programs and through the event reconstruction 
software. 

The radiative corrections were calculated to be (10-15)% for the selection cuts used in the 
analysis and for the Q 2 and W dependences found in the data. They are taken into account in 
the cross sections given below. 

5 Analysis and cross sections 
5.1 Data selection 

For the selection of exclusive p° candidates, the off-line analysis required: 

• a scattered electron energy, as measured in the calorimeter, greater than 5 GeV. The elec- 
tron identification algorithms used in this analysis were optimised to have high efficiency 
(>97%); 

• 5 = J2i Ei(l — cos#i) > 35 GeV, where the sum runs over all calorimeter cells; this cut 
reduces the radiative corrections and photoproduction background; 

• two tracks with opposite charge, both associated with the reconstructed vertex; if there 
was a third track at the vertex, it should be from the scattered electron. Each of the 
two tracks was required to have a transverse momentum above 0.16 GeV and a polar 
angle between 25° and 155°; this corresponds to the region where the CTD response and 
systematics are well understood; 

• a measured vertex (Z vtx ), as reconstructed from VXD and CTD tracks, to be in the range 
-50 < Z vtx < 40 cm; 

• events with a scattered electron whose impact point in the RCAL was outside the square 
of 32 x 32 cm 2 centered on the beam axis or events with an RHES impact point outside 
the square of 26 x 26 cm 2 ; this requirement controls the determination of the electron 
scattering angle; and 

• the residual calorimeter energy not associated with the electron to be compatible with 
the p° momentum measured in the tracking system, E p CAL jP p < 1.5 (see Fig. |l]a), where 
E P CAL is the calorimeter energy excluding that due to the scattered electron and P p is the 
sum of the absolute values of the momenta of the two oppositely charged tracks. This 
cut suppresses backgrounds with additional calorimeter energy unmatched to the tracks 
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and events with proton dissociation depositing energy in the calorimeter towers around 
the FCAL beampipe. Also shown in Fig. |l|a is the distribution from the p° DIPSI Monte 
Carlo events, indicating that only a small fraction of the exclusive p° events are removed 
by this cut. Fig. [l]b shows the same distribution after the final selection indicating good 
agreement with the expected distribution. 

A total of 352 events passed these selection requirements. 

Possible backgrounds to the exclusive reaction (1) are from p° events with additional unde- 
tected particles, from and u production and from proton dissociation events where Mx is 
small and therefore does not deposit energy in the calorimeter. To reduce these backgrounds, 
two additional cuts were imposed: 

• 0.6 < M^+t,.- < 1.0 GeV; this selection reduces the contamination from and u produc- 
tion in the low 7r + 7r~ mass region as well as higher mass resonant states and non-exclusive 
events in the high mass region; and 

• Pt < 0.6 GeV 2 ; this cut reduces non-exclusive background and proton dissociation events. 
Fig. [I|c shows the measured, uncorrected p\ distribution for the selected p° events, indi- 
cating a clear excess of events above a single exponential for p\ > 0.6 GeV 2 , consistent 
with the presence of proton dissociation events which, in hadron-hadron scattering, have 
a less steep p\ distribution pjfl . (The acceptance is relatively flat, rising by about 5% 
from p 2 T = to 0.6 GeV 2 .) 

Fig. |l]d shows a scatter plot of Q 2 versus x for the 140 events which pass the above criteria. 
The efficiency drops sharply at small Q 2 (due to the cut on the electron impact point in 
RCAL) and at small and large y (due to the requirements on the 7r + ,7r~ tracks). To remove 
poorly reconstructed events and to select a region of phase space where the acceptance is well 
determined and relatively constant as a function of the kinematic variables, two additional 
kinematic cuts were applied to the data: 

• 7 < Q 2 < 25 GeV 2 and 0.02 < y < 0.20. 
The final p° sample contains 82 events. 



5.2 Background estimates and acceptance corrections 

The DIPSI p° Monte Carlo simulated events were used to correct the data for acceptance and 
detector resolution. The acceptance (which includes the geometric acceptance, reconstruction 
efficiencies, detector efficiency and resolution, corrections for the offline analysis cuts and a 
correction for the M n + n - cut) in this region of Q 2 varies between 40% and 55%. The acceptance 
is constant at about 47% as a function of y, p\ or M n + n - in the above kinematic region. The 
resolutions in the measured kinematic variables, as determined from the Monte Carlo events, 
are 6% for Q 2 and 2% for y. 

Fig. |^a, which shows the uncorrected 7r + 7r~ mass distribution for the events passing all of 
the final cuts (except for the M n +^- cut, but with a M K + K - > 1.05 GeV cut, when the tracks 
are assigned a kaon mass, to remove <fi events), indicates a pure sample of p° events. A p° 
non-relativistic Breit-Wigner form, with a constant background, is fit to the mass spectrum 
between 0.6 and 2.0 GeV. The resulting parameters for the p° mass and width are 774 ±9 MeV 
and 134 ± 20 MeV, respectively, to be compared with the values 769.9 MeV and 151.2 MeV 
from the Particle Data Group [f2{|. The fit also includes a flat background estimate of (4 ±4)% 
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for M^+^- masses between 0.6 and 1.0 GeV, as determined from comparing the Monte Carlo 
p° events with the data for M n + n - masses between 1.0 and 1.5 GeV. Background contributions 
from exclusive <fi and u events were estimated to be at the 1% level and are included in the 
above background estimate. 

Since the proton was not detected, the proton dissociation background contribution had to 
be subtracted. This was done using the RHODI event generator combined with the detector 
simulation. The normalisation was obtained by requiring that the Monte Carlo generated 
sample have the same number of events with energy between 1 and 20 GeV in the FCAL 
as for the data (7 events) when the constraint that E p CAL jP p < 1.5 was relaxed to [E P CAL — 
Efcal) I ' P P < 1-5 and the additional constraint 9 n ± > 50° was imposed. Assuming an M x 
dependence of the form l/M 2 ; 25 , as measured by the CDF experiment [30] for pp — > p + X, 
yielded a contribution of (22 ± 8 ± 15)% where the systematic error was obtained from varying 
the exponent of 1/Mx between 2 and 3. This is consistent with an estimate from the excess 
above the exponential in the p\ distribution mentioned above. In the exclusive p° sample 
under study here, there are no events with an electron energy in the range 5 < E' e < 14 GeV 
and so the photoproduction background is negligible. No events were found from the unpaired 
bunches demonstrating that the beam-gas background is also negligible. The overall background 
contamination was estimated to be A = (26 ± 18)%. Unless explicitly stated otherwise, this 
background was subtracted as a constant fraction for the cross sections given below. 



5.3 The ep cross section 

The cross section, measured in the kinematic region defined above, is obtained from a(ep — > 
ep°p) = N(l — A)(7i/(C/2 • A ■ L int ), where N (= 82) is the observed number of events after 
all cuts with 0.6 < M nn < 1.0 GeV, A is the background estimation, A is the acceptance as 
discussed above, L int is the integrated luminosity of 0.55 pb -1 and C<i is the correction for QED 
radiative effects. These radiative corrections were calculated for the exclusive reaction using 
the x and Q 2 dependences found in this experiment (see section 6.1) and vary between 1.10 (at 
low Q 2 and low y) and 1.14 (at high Q 2 and high y). A systematic uncertainty of ±0.10 was 
included on this correction to account for the uncertainties in the cross section dependences 
on x and Q 2 . To compare later with results from the NMC experiment, which has determined 
exclusive p° cross sections integrated over all p\ 0, C\ is a 4.5% correction for the cross section 
in the p\ range between 0.6 and 1.0 GeV 2 based on the slope of the distribution measured in 
the present analysis. The corrected ep cross section for exclusive p° production at yfs = 296 
GeV is 

a(ep -> ep°p) = 0.21 ± 0.03(stat.) ± 0.06(syst.) nb, 

integrated over the ranges 7 < Q 2 < 25 GeV 2 , 0.02 < y < 0.20 and p 2 T < 1.0 GeV 2 , with 
acceptance corrected < Q 2 > and < W > of 11.0 GeV 2 and 78.9 GeV, respectively. 

The quoted systematic uncertainty is derived from the following (the systematic uncertainty 
for each item is indicated in parentheses): 

• the cuts used to remove non-exclusive backgrounds were varied and independent analyses 
using differing selection cuts and background estimates were compared to the previously 
described analysis: tracks were matched to the calorimeter energy deposits and events 
containing an energy in excess of that of the p° of more than 0.4, 1 or 2 GeV were 
discarded; events were selected based on the position of the electron measured by the 
calorimeter rather than by the RHES and the cut on the impact position of the electron 
in RCAL was varied (10%); 
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• using different trigger configurations (8%), 



• the cuts on the tracks were varied. The lower cut on the transverse momentum was varied 
between 0.1 and 0.2 GeV and different polar angle selections were made. The maximum 
variation occurred for p T > 0.2 GeV (9%); and 

• events from different Monte Carlo generators were used to calculate the acceptance and 
efficiency (7%). 

Adding these in quadrature to those from the uncertainties due to background subtraction 
(24%), luminosity (3.3%) and radiative corrections (10%) yields 31% as the overall systematic 
uncertainty. 



5.4 The 7*j9 cross sections 

The ep cross section was converted to a 7*p cross section as follows. The differential ep cross 
section for one photon exchange can be expressed in terms of the transverse and longitudinal 
virtual photoproduction cross sections (see @) as: 



d 2 a(ep) a 



dxdQ 2 2ttxQ 2 



i + (i - yf) ■ 4 p {y, Q 2 ) + 2(1 - y) ■ aT'iv, Q' 2 



The virtual photon-proton cross section can then be written in terms of the electron-proton 
differential cross section: 

a(lP^p P ) = (a T +ea L ) = - — 2 (7) 

where Ft, the flux of transverse virtual photons, and e, the ratio of the longitudinal to transverse 
virtual photon flux, are given by 

«(1 + (!-»)') md 2(1-,) 



2ir.iQ 2 (1 + (1 - yf) 

Throughout the kinematic range studied here, e is in the range 0.97 < e < 1.0. 

Using Eq. (0), cr(-y*p — > p°p) was determined with the flux calculated from the Q 2 , x and y 
values on an event-by-event basis. The 31% overall systematic uncertainty on a(ep) applies to 
every value for a(^*p — > p°p) and thus becomes an overall normalisation uncertainty. 



6 Results 

6.1 Q 2 and p\ distributions 

After correcting for detector acceptance and backgrounds, the cross sections were obtained 
as a function of Q 2 . Fig. [|b displays the Q 2 dependence of the 7*p — > p°p cross section for 
events in the x range between 0.0014 and 0.004. Also displayed in Fig. |2|b are data from the 
NMC experiment . The ZEUS values of the cross sections are larger than those of the NMC 
experiment. However, it should be noted that for this figure as well as for Figs. 2d, 3 and 4 
the different experiments have different mixtures of longitudinal and transverse photon fluxes 
(e varies from 0.5-0.8 for the NMC results). More importantly, the region of 7*p cm. energy 
of the NMC experiment (8-19 GeV) is lower than that in this experiment (40-130 GeV). 



8 



A fit of the form Q~ 2a to the distribution of the ZEUS data in Fig. |2|b yields the power 
of the Q 2 dependence. To study the systematic uncertainty on the value of a, a maximum 
likelihood analysis was performed with the cross section factorised as: 

v(l*P^P°p) ~ {Q 2 T a ■ x~P ■ e~ bp T . (8) 

This study, applied to the 82 events in the final data sample in the region 0.02 < y < 0.20, 
p\ < 0.6 GeV 2 and 7 < Q 2 < 25 GeV 2 , yields results similar to those obtained from fitting 
Fig. ^|b. The best estimate of the Q 2 dependence is 2a = 4.2 ± 0.8(stat.)1lQ^(syst.), where 
the systematic uncertainty comes from the variation in the value of a obtained from the two 
different fitting methods as well as from the variation obtained from the systematic studies 
described in section 5.3. 

The uncorrected p\ distribution was presented in Fig. and showed an exponentially falling 
behaviour, with an excess of events for p\ > 0.6 GeV 2 , as discussed previously. After correcting 
for detector acceptance and resolution, a fit in the range < p\ < 1.0 GeV 2 of the form: 

<h Ael-*&> + Be<~te\ (9) 



2 



dp' 

was performed. In Eq. (|9|) the contribution from the second term, which was constrained to be 
22% of the number of events for p\ < 0.6 GeV 2 , represents the proton dissociative background 



contribution which was assumed to have a slope half that of the exclusive reaction |pq| . The fit 
yields a slope parameter of b = 5.1^Q^(stat.) ± 1.0(syst.) GeV" 2 , which is consistent with that 
found in the maximum likelihood fit. The systematic uncertainty comes from the variation in 
the value of b obtained from fits without the second term in Eq. (|9|), the maximum likelihood 
fit and from the systematic studies described in section 5.3. This value of b is about half that 
found in elastic p° photoproduction [|1(J and is in agreement with the result from the NMC 
experiment 0. 

6.2 p° decay distribution 

The p° s-channel helicity decay angular distribution H(cos9h, (fih, &h) can be used to determine 
the p° spin state fl3l| , where 9h and (fih are the polar and azimuthal angles, respectively, of the 



7r + in the p° cm. system and $/,, is the azimuthal angle of the p° production plane with respect 
to the electron scattering plane. The quantisation axis is defined as the p° direction in the 7*p 
cm. system. Only the cos^ dependence is presented here. After integrating over <ph and 
the decay angular distribution can be written as: 

1 f/ * V V^ + (3^-l)cosX], (10) 



N dcos9 h 4 1 

where the density matrix element Tqq represents the probability that the p° was produced 
longitudinally polarised by either transversely or longitudinally polarised virtual photons. 

The helicity cos9h distribution (uncorrected for background, since the dominant contribution 
to the background is from proton dissociation which is expected to have the same helicity as 
the p°p final state) is shown in Fig. [|c. The curve shown in the figure is from a maximum 
likelihood fit to the form of Eq. (0) yielding = 0.6 ± 0.1±g;f at < Q 2 > = 11.0 GeV 2 and 
< W > = 78.9 GeV, where the first uncertainty is statistical, and the second is derived from 
the variations of the result when different ranges in cos#h were used in the fit and when the 



9 



systematic studies of section 5.3 were used. In Fig. ^|d this measurement of is compared to 
other published data at various values of Q 2 . These data show the presence of both transversely 
and longitudinally polarised p°'s at high Q 2 (above 2 GeV 2 ). If SCHC is assumed, an estimate 
of R, the ratio of longitudinal to transverse cross sections, for p° production is obtained 0: 

P — ° L — 1 r 00 _ 1 c+2.8 

cjt e i — r 00 

(where the statistical and systematic uncertainties in t[]q have been added in quadrature). This 
may be compared with the value of R = 2.0 ± 0.3 at < Q 2 > = 6 GeV 2 and < W >~ 14 GeV 
from the NMC experiment 0]. 

6.3 The W and x dependences of the 7*p — > p°p cross section 

Fig. |^ shows a compilation [2,4,5,7-9] of photoproduction and selected leptoproduction exclusive 
p° cross sections as a function of both Q 2 and W. In this figure the cross sections obtained 
in this analysis are shown as a function of W at Q 2 = 8.8 and 16.9 GeV 2 . The cross sections 
at different W values (and slightly different < Q 2 >) were scaled to Q 2 = 8.8 and 16.9 GeV 2 
using the measured Q~ A 2 dependence in order to compare with the NMC cross sections^ from 
deuterium at the same values of Q 2 . At high energies, W > 50 GeV, data exist only at Q 2 = 
0, 8.8 and 16.9 GeV 2 . The real (Q 2 = 0) 'jp 'elastic' cross section [10| shows only a slow rise, 
consistent with that seen in the photon-proton total cross section. At small Q 2 (< 2.6 GeV 2 ), 
the data first decrease with increasing W followed by a slow increase. No high energy data are 
yet available to see how the increase develops. At higher Q 2 , the cross sections rise strongly 
with increasing W. At Q 2 = 8.8 GeV 2 and W ~ 100 GeV, the cross section is about a factor 
of six larger than at W = 12.9 GeV [|7J. This strong increase in the 7*p — > p°p cross section is 



in contrast to that expected from the Donnachie and Landshoff model |TT], |12| . 

To compare with the QCD calculations of Brodsky et al. and Donnachie and Landshoff, the 
ZEUS and NMC cross sections are shown as a function of x at Q 2 = 8.8 and 16.9 GeV 2 in Figs. 
|a, b. The total cross sections are predicted from the pQCD calculations of Brodsky et al. 
using the longitudinal contribution to the differential cross section at t — 0, (see Eq. 



da(x,Q 2 ,p 2 T ) ^ 2 _ (1-e b ) d(a T + a L ) 
a{x, ) - I ap T 



[} ~ e- h ) 1 | da, 



R ' dt 



dp 2 T b dt t=Q ., „ t=0 

(11) 

where b is the slope of the p\ distribution. The measured values of e = 0.98, R = 1.5 and 
b = 5.1 GeV~ 2 were used to calculate a(x, Q 2 ). For the gluon momentum density, xg(x,Q 2 ), 
the form xg ~ x s ■ (1 — x)' n was assumed. The values of 5 and i] were allowed to vary within 
the ranges (—0.25 to —0.39) and (6.44 to 4.81) respectively. These ranges were determined 
from the leading order gluon density extracted from the scaling violations of the ZEUS 
F 2 measurements [jlTJ. The uncertainty in Eq. (|TT|) arising from the la range in the gluon 
density |18[] is shown as the light shaded region in Fig. ^. The uncertainty in the prediction 



arising from measurement uncertainties of R and b when added in quadrature to that of the 
gluon distribution yields the larger dark shaded area. The shaded areas in Fig. f| are restricted 
to x < 0.01 where Eq. ([3D is valid. The range in the predicted cross sections at a given x 
is dominated by the uncertainty on 5. Since the calculation is made in DLLA, the value of 

4 Since the EMC and NMC data cover approximately the same kinematic region, the more recent NMC 
data|Q have been chosen to make comparisons. 
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Q 2 at which the gluon density and a s are evaluated is not defined to better than a factor of 
two. Varying the Q 2 scale for a s (Q 2 ) ■ xg(x,Q 2 ) from Q 2 /2 to 2Q 2 changes the prediction by 
about 50%. The x scale in the calculation can also range from x/2 to 2x so that the curves 
can be shifted left or right to reflect this uncertainty. At the present level of precision of the 
measurement and theory, the data are consistent with the pQCD calculation of Brodsky et al. 
The hatched region shows the range of the predictions of Donnachie and Landshoff based on 
soft pomeron exchange [II] with o = | ^ ^ . The range of the hatched area comes from the 
uncertainty in the measured value of b. The data do not agree with these expectations, being 
typically a factor of three above the predictions. 



7 Conclusions 

Exclusive p° production has been studied in deep inelastic electron-proton scattering at large Q 2 
(7 - 25 GeV 2 ) in the 7*p centre of mass energy (W) range from 40 to 130 GeV. Cross sections 
are given for both the ep and 7*p processes. The cross section for the 7*p process exhibits 
a Q~( 4 - 2±a8 -o.5) dependence. The *y*p — > p°p cross section at these large Q 2 values shows a 
strong increase with W at HERA energies over the lower energy NMC data, in contrast to the 
behaviour of the elastic photoproduction cross section. Both longitudinally and transversely 
polarised p°'s are produced. The Q 2 dependence, the polarisation and the slope of the p\ 
distribution are consistent with those observed at lower energies. However, the cross sections 
are significantly larger. The Donnachie and Landshoff prediction for soft pomeron exchange 
underestimates the measured cross sections while the data are consistent with the perturbative 
QCD calculation of Brodsky et al. given the present knowledge of the gluon momentum density 
in the proton. 
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Figure 1: (a) The distribution of E p CAL jP p for the candidate events; (b) the same distribution 
for the final sample of 82 events; (c) the p\ distribution; and (d) a scatter plot of Q 2 versus 
x for the selected p° events. These plots are not corrected for detector and trigger efficiencies 
and acceptances. The histograms in (a-c) are obtained from the DIPSI p° Monte Carlo sample 
after detector and trigger simulation. In (d) the lines correspond to the region in Q 2 and y 
selected for this analysis. 
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Figure 2: (a) The 7r + 7r~ invariant mass distribution for the final sample of events; the curve 
is a maximum likelihood fit with a non-relativistic Breit-Wigner distribution plus a flat (4%) 
background (see text for details); (b) the cross section for 7*p — > p°p as a function of Q 2 for 
0.0014 < rr < 0.004. Also shown are data from the NMC experiment [7]; the errors shown are 
just the statistical errors. The ZEUS (NMC) data have an additional 31% (20%) normalisation 
uncertainty (not shown); (c) the cos^^ distribution for the decay 7r + , in the s-channel helicity 
system, corrected for acceptance, for 7r + 7r~ pairs in the mass range 0.6-1.0 GeV. The curve is a 
fit to the form of Eq. (10); (d) the p° density matrix element, Tqq, compared with results from 
fixed target experiments [2,3,7] as a function of Q 2 . The thick error is the statistical error and 
the thin error is the systematic error added in quadrature. 
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Figure 3: The YP ~^ P°P cross section as a function of W, the 7*p centre of mass energy, for 
several values of Q 2 . The low energy data (W < 20 GeV) come from fixed target experiments 
[2,4,5,7,8]. The high energy data (W > 50 GeV) come from the ZEUS experiment [10] and the 
present analysis. The ZEUS data at Q 2 = 8.8 and 16.9 GeV 2 have an additional 31% systematic 
normalisation uncertainty (not shown); the data from Refs. [2], [4] and [7] have additional 10%, 
25%, and 20% normalisation uncertainties, respectively. 
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Figure 4: (a) The cross section, a{^*V ~~ > P°P)-> as a function of x at a value of Q 2 = 8.8 
GeV 2 . (b) A similar plot for data at Q 2 = 16.9 GeV 2 . The errors shown are only the statistical 
errors. In addition, there is a 31% systematic uncertainty which is not shown but applies 
to the overall normalisation. Also shown is the NMC result (which has an additional 20% 
normalisation uncertainty [7]). The shaded area corresponds to the predictions of Eqs. (3) 
and (11) for x < 0.01. The range of the predictions shown by the light shaded area is a result 
of the experimental uncertainty of the gluon distribution [18]. The dark shaded area includes 
the uncertainties on b and R added in quadrature to that of the gluon. In addition, there is a 
50% uncertainty in the predicted cross section from the choice of the Q 2 scale; furthermore, the 
value of x in Eq. (3) is only defined to within a factor of 2. The hatched area displays the cross 
section expected from the soft pomeron model [11]. The range comes from the uncertainty in 
the measured value of b. 
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